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摘 要 :大 多 数 悬 臂 梁 的 研究 都 是 基于 确定 性 振动 分 析 或 将 外 部 激励 假定 为 高 斯 过 程 , 鲜 少 涉及 非 
高 斯 激励 。 针 对 非 高 斯 激励 下 带 附加 质量 的 悬臂 梁 的 响应 特征 进行 研究 。 首 先 ,通过 Lagrange 方 
程 建立 考虑 曲率 非 线 性 的 悬臂 梁 振 动 方程 ;然后 ,通过 指数 多 项 式 闭合 法 求解 结构 响应 的 概率 密度 
函数 ,并 分 析 了 平均 到 达 速 率 、 激 励 强度 及 附加 质量 对 结构 响应 的 影响 ;同时 分 析 了 非 平 稳 非 高 斯 
激励 下 非 线性 结构 参数 对 响应 的 影响 。 计 算 结果 表明 ,指数 多 项 式 闭合 法 与 模拟 解 吻 合 很 好 ,平均 
到 达 速 率 越 大 ,结构 响应 的 概率 密度 函数 曲线 收缩 越 明 显 。 
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Abstract: Most of the research on cantilever beam is based on deterministic analysis or the assumed 
Gaussian excitation, non-Gaussian excitation is rarely involved. Stochastic responses of cantilever beam 
with lumped mass subjected to non-Caussian random excitation are studied in this paper. Motion equation 
of the structure ls established by the Lagrange’s equation which considers curvature nonlinearity , approxi- 
mate transient solutions of responses of structure are obtained by exponential polynomial closure (上 PC ) 
method. Effects of average arrival rate ,excitation intensity and lumped mass on the structural response are 
analyzed. Meanwhile ,the influence of nonlinear structural parameters on structural response under non-sta- 
tionary and non-Caussian excitation is also analyzed. Numerical results show that results obtained by the 
EPC method agree well with simulation solution. With the increase of the average arrival rate ,the probabi- 


lity density function curve of structural response presents centralization. 
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很 多 结构 可 以 简化 成 带 附加 质量 的 悬臂 深 模 
型 ,如 高 答 建 筑 海上 钻井 平台 、 信 号 塔 和 交通 灯 框 
架 等 ,目前 已 有 一 些 针对 此 类 模型 的 研究 。CRESPO 
DA SILVA 通过 广义 Hamilton 准则 推导 了 无 轴 向 变 
形 悬 臂 梁 的 三 次 非 线性 动力 学 方程 ,并 考虑 了 曲 
率 非 线性 ?1 。 随 后 ,HYFER 研究 了 悬臂 梁 在 受 迫 振 
动 下 的 空间 动力 学 行为 5] 。 基 于 Timoshenko 梁 准 
则 ,HAMDAN 建立 了 带 附 加 集中 质量 块 的 基 臂 粱 自 
由 振动 方程 个。CHENG 利用 Hamilton 准则 和 La- 
grange 方程 推导 了 考虑 系统 阻尼 与 高 斯 白 品 声 激励 
的 上 甚 臂 梁 运动 方程 ,并 利用 随机 等 价 线 性 化 法 研究 
了 结构 响应 的 统计 和 矩 5 。WANG 推导 了 旋转 锥 形 
悬臂 梁 的 横向 振动 方程 “”。 目 前 ,大 多 数 的 研究 是 
基于 确定 性 激励 或 者 高 斯 激励 ,而 基于 非 高 斯 激励 
下 芒 壁 梁 的 响应 特征 却 鲜 有 研究 。 众 所 周知 ,如 果 
系统 响应 是 马尔 科 夫 过 程 , 则 系统 响应 的 概率 密度 
国 数 满足 广义 Fokker-Planck-Kolmogorov (FPK ) 方 
程 。 而 非 高 斯 激励 下 系统 的 FPK 方程 含有 对 响应 
的 四 阶 偏 微分 项 ,这 使 得 它 比 高 斯 激励 下 非 线性 系 
统 的 啊 应 求解 过 程 要 复杂 得 多 ,计算 效率 显著 降低 。 
然而 众所周知 ,自然 界 中 大 多 数 激 励 的 脉冲 幅 值 都 
不 是 连续 的 ,如 罕 遇 地 震 的 地 面 运动 加 速度 、 冲 击 
波 、 作 用 在 离 岸 结构 海上 的 波浪 、 路 面 或 轨道 的 不 平 
坦 对 车 辆 的 作用 等 ,这 些 激励 相 邻 的 脉冲 幅 值 间 具 
有 一 定 的 时 间 间 隔 , 简单 地 将 外 部 激励 看 作 高 斯 白 
噪声 激励 将 使 得 计算 的 响应 概率 密度 子 数 偏 于 不 安 
全 。 因 此 ,很 多 国内 外 学 者 研究 了 非 高 斯 激励 下 系 
统 响应 的 求解 方法 。 

目前 很 多 学 者 提出 了 非 高 斯 激励 下 系统 反应 分 
析 的 方法 , 如 随机 等 价 线性 化 法 、 随 机 平均 
法 "Wiener 路 径 积分 法 "及 指数 多 项 式 闭合 法 
(EPC) 等。 这些 方法 大 多 是 从 高 斯 激励 下 系统 
反应 分 析 方 法 发 展 而 来 。 在 众多 方法 中 ,指数 多 项 
式 闭合 法 具有 和 较 高 的 精度 与 计算 效率 ,而 且 可 以 得 
到 系统 瞬 态 响应 的 近似 解析 式 。 因 此 ,本 研究 首先 
通过 Lagrange 方程 建立 非 高 斯 激励 下 带 附加 集中 质 
量 甚 辟 梁 的 运动 微分 方程 ,利用 指数 多 项 式 闭 合法 
求解 结构 响应 的 概率 密度 函数 ,并 分 析 非 高 斯 激励 
的 平均 到 达 速 率 和 激励 强度 对 结构 响应 的 影响 、 附 
加 质量 对 结构 参数 及 响应 的 影响 。 同 时 ,分 析 了 非 


平稳 非 高 斯 激励 下 激励 参数 对 响应 的 影响 。 


1 随机 微分 方程 


带 附 加 质量 的 悬臂 梁 如 图 1 所 示 。 质 量 块 履 距 
基础 4 的 距离 记 为 ,mm 为 杆 件 的 线 密 度 ,假定 杆 件 
是 均匀 的 , 则 观 是 一 个 常数 , 是 质量 块 必 的 转动 惯 
量 , a,(1) 是 地 面 运动 加 速度 。 


| CAG) 
Em 


图 1 悬臂 粱 示意 图 
Fig.1 Structure sketch map under external excitation 
通过 Lagrange 方程 可 以 建立 悬臂 粱 的 运动 方 
程 。 采 用 文献 [5 ] 的 结果 ,系统 的 势能 可 以 表示 为 
U = Fen | x(aé = NE + (I 
(1) 


其 中 , x 表示 曲率 。 根 据 Euler 梁 理 论 ,考虑 附加 质 
量 的 影响 ,该 系统 的 动能 可 以 表示 为 


大半 型 | (i 1 7) 和 E+ 工 永 
2 儿 2 
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其 中 :了 Et) = y(é,t) ZE = 2(é,t) - a(t) 。 

因此 ,Lagrange 方程 可 以 表示 为 
Po re (55 
该 方程 里 含有 响应 与 时 间 的 偏 导 项 ,很 难 求 得 
系统 响应 的 解 。 通 常 将 响应 分 解 成 状态 空间 与 时 间 

相互 独立 的 量 , 即 

y(E,t) = P(E)x(L) (4) 
其 中 : p(E) 表示 梁 的 一 阶 模 态 函数 ; x(1) 为 模 态 坐 
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标 。 一 般 情 况 下 ,悬臂 粱 的 低 阶 模 态 起 决定 性 作用 ， 
因此 ,本 研究 仅 考 虑 一 阶 模 态 下 悬臂 粱 的 振动 方程 。 

为 方便 起 见 ,悬臂 粱 的 附加 集中 质量 块 设置 在 
深 的 末端 。 根 据 文 献 [14 ] , 模 态 郧 数 可 以 取 
oD(E) = coshyé - cosvé -CsinhyE - sinvé)] 
(5 ) 
对 于 悬臂 梁 , 取 > =0.734 095 5,Z =1.875 104 1。 
模 态 坐标 下 的 Lagrange 算 子 可 以 表示 为 
= Di -2d + Si -Bao 
6 a 
Ve 
其 中 :6 = EL/(ml);a = J/(me);a, = M/(m)); 
CQ1 ~Qs 表示 为 


x + (a, + Do (6) 


1 
oy = | pd t ag D+ aap Cn) 


1 
Q2 一 | ea 
1 € 2 
a fe ee 
a (7) 
a 和 “de 
: 有 é=n 
1 
ay = | op" 上 
1 
Qs = | ed + Qa.p (nm) 
将 式 (6) 代 入 Lagrange 方程 ,得 
a +6 ox +20° Be + 人 (4 i 
Qs a,(t) (8) 


考虑 系统 阻尼 ,悬臂 粱 的 运动 方程 可 表示 为 


X+HX+3Y+oac +BOxi2+i22) = se av(i) 


(9) 
其 中 

a a 28’ 

三 < 2 2 -2% 4 0 _ 

p= 0 = = ga = 区 全 ， 

Q3 Qs 

= 一 10 
B ee (10) 


准确 描述 外 部 激励 对 随机 振动 分 析 具 有 重要 意 
义 。 多数 文献 表明 地 震 激 励 a,(t) 是 离散 不 规则 过 
程 … ,每 次 观测 得 到 的 时 间 历 程 样本 都 不 相同 ， 
但 总 体 样本 符合 泊 松 过 程 ,而 激励 样本 的 功率 谱 密 
度 大 致 是 均匀 分 布 的 ,因此 可 以 用 非 平稳 泊 松 白 噪 
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声 模型 描述 , 即 平稳 泊 松 白 品 声 激励 WW(1) 与 时 间 
调制 函数 h(t) 相 乘 。 假 定 V 为 泊 松 计数 过 程 ,其 到 
达 速 率 函数 A(1) 通常 取 一 段 时 间 内 脉冲 振幅 的 平 
均 达 到 速率 A 。7, 是 在 7; 时 刻 一 系列 零 均 值 独立 同 
分 布 的 脉冲 幅 值 ,因此 , WW(1) 可 以 表示 为 


N 
W(t) = 》76G -7) (11) 
$= 
其 中 , 6( .) 是 Dirac delta 函数 。 


考虑 式 (9) 中 含有 位 移 与 加 速度 耦合 项 ,因此 , 需 
要 采用 近似 方法 将 振动 方程 处 理 。 式 (9) 可 整理 为 


1 . eh(t) W(t) 
三 3 2)+ 一 -一 

1 +Bx MX+YyX+oAax +BxxX 1 + Br 
(12) 


假定 Bx” < 1 ,1A(1 + Bo) 可 以 取 其 泰勒 级 数 
展开 式 的 前 两 项 , 即 


1 2 
TB = 1 — Bx (13) 
因此 , 式 (12) 可 以 近似 写 为 
YY = 一 (1 — Bx ) Cx + yx 十 ao + Bxx ) 十 
g(1 _ BoACDWVC) (14) 


令 x1 = x,x， = ,可 以 将 式 (14) 写 成 如 下 的 状 


Wi 
x, = 一 (1 — Bx1) Cx + YX + Oxi + Bxix2) + 
eGCx h(t) W(t) (15) 


其 中 , G(x1) = ae(1 - Bxi) ,初始 条 件 为 
E[xi;0] = EL[x,;0] = 0,E[%?;0] = o?(0), 
E[xi;0] = o3(0) 


2 指数 多 项 式 闭 合法 


在 泊 松 白 噪 声 激 励 下 ,悬臂 梁 的 响应 是 马尔 科 
夫 过 程 ,其 概率 密度 函数 满足 如 下 FPK 方程 


3 0 
+x, 一 一 (1 -Bxi) (x2 + YX + Ox] +Brit2) 一 
1 2 


0[(1 -Bxi) (yx, +yai +axl +Bxix2)] 


0xX2 


1 07 
7ABLY IG (Cx) (+ 
二 2 

1 9 
1AELY GC (x ) h(t) 
5 OX? 


4 
ye =0 (16) 
41! OX? 
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其 中 ,p 是 响应 xi wx, 的 联合 概率 密度 函数 。 通 常 方 
程 仅 保留 非 高 斯 激励 的 前 四 阶 项 ,这 是 因为 高 阶 项 
对 整个 方程 的 贡献 很 小 。 

该 方程 是 一 个 二 维 演 化 类 非 线性 偏 微分 方程 ， 
它 的 解析 解 很 难 求 出 。 通 常 采用 数值 方法 近似 计算 
系统 响应 的 概率 密度 函数 。 假 定 近似 解 表达 式 为 
plx;a(i)]= Cexp{Q,[x;a(t)]} (17) 
其 中 , C 是 归 一 化 常数 ,而 
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N, 
QO,[x;a(i)] = > a (Fr tr 
f= 


oe i 
rkl+r [kl= 1,2,.…,n 
将 式 (17) ~ (18) 代 入 式 (16) ,等 式 右 侧 则 出 现 
误差 , 即 


30, 00, 3 90, oO[ (1 -Bxi) (pxs +Bxixz)] 
Alxsa()1 = 全 + 和 (Be ty to 45- 2 人 
人 2 .100, 1/90, 1 ， Po A A A 
T7AELY GC (0)h (| ee | | ] -21c | 0 
0.0, 00 2 0 0, 190 311- _ 
3 | | | | ]}itx;al)] =R[x;a(t) ]b[x;a(t)] (19) 
采用 人 徊 辽 金 法 优化 RLx;a(1)] 使 误差 最 小 , 即 [= [Ca (22) 
0 


sRixsa() lo,lx;o(t) ldr = 0， 
r = 1,2,.,N, (20) 
其 中 w,[x;og(1)1 是 权重 函数 ,可 以 取 
em =2 wf [x ,ot) fl ,0(t) ] 


ml[r],m,[r] =0,1,2,.…,n 


ml[rj+m[r]=1,2,.…,n, r=1,2,.…,N 


(21) 
其 中 : fi[xj ,oj (2)1, 记 Lx ,0;(1)] 分 别 是 由 高 斯 闭 
合法 得 到 的 x 和 x%, 的 边缘 概率 密度 函数 ; o1 (7?)、 
0,(t) 分 别 是 x 和 2 在 时 刻 1 的 方差 。 

经 过 整理 , 式 (20) 变 成 了 以 待定 系数 w 为 变量 的 
常 微分 方程 组 ,通过 四 阶 龙 格 库 塔 微 分 程序 ,可 以 较为 
快速 地 得 到 待定 系数 的 解 。 求 出 待定 系数 后 ,通过 归 
一 化 条 件 , 可 以 得 到 系统 响应 的 联合 概率 密度 函数 。 


3 算 例 分 析 与 讨论 


悬臂 梁 采 用 长 度 / =6 m 的 B32a( GB/T 706 - 
1988 ) ,计算 方向 为 截面 弱 轴 ,一 阶 模 态 阻尼 比 = 
0. 02 。 考 虑 两 种 不 同 的 情况 , a, = o = 0 及 a = 
o 天 0 。 前 者 可 以 模拟 高 管 建筑 .信号 塔 海上 钻井 
平台 等 ,后 者 可 以 模拟 交通 灯 框 架 广告 牌 等 。 本 研 
究 采 用 的 单位 :长 度 单位 取 x 10 m 和 时 间 单 位 取 
x0.1s, 算 例 的 结果 均 是 悬臂 梁 在 端点 处 的 响应 统 
计量 。 算 例 考 虑 了 平稳 泊 松 白 噪 声 和 强度 非 平稳 泊 
松 白 噪声 激励 的 情况 ,考虑 一 致 性 ,引入 Arias 强 
度 , 即 


衡量 地 面 运 动 的 强度 。 同 时 ,为 了 验证 该 计算 
方法 的 精度 ,采用 蒙特 卡 洛 模拟 解 与 该 近似 解 作 对 
比 ,样本 数量 取 10”。 由 于 抗震 分 析 一 般 基于 承载 力 
与 变形 计算 ,因此 ,本 研究 仅 展 示 结 构 的 位 移 响应 结 
果 。 对 于 速度 响应 的 结果 ,限于 篇 幅 , 在 此 不 予 展示 。 


3.1 Box-Car 型 


首先 考虑 平稳 泊 松 白 噪 声 激 励 下 的 情况 , 即 
h(t) = 1 (23 ) 
取 到 达 率 A = 20 ,激励 强度 为 AELY] = 
0.5 m*/s'。 不 考虑 附加 质量 , 即 a = a,，= 0 。 通 过 
高 斯 闭合 法 ,可 得 到 响应 的 统计 矩 ,结果 如 网 2 所 
示 。 其 中 ,响应 初 值 取 oC(0) = 201(0) = 0.001 。 


Sa El[x’ (EPC) Elxi](MCS) 
“Elxx,](EPC) ©——E[Lxx,](MCS) 
0.006 -—-—-ELx:](EPC) E[Lx;](MCS) 


20 40 60 80 100 
1(10 s) 


图 2 平稳 激励 下 结构 响应 的 统计 和 矩 
Fig.2 Moments of system responses under stationary excitation 
从 图 2 可 知 ,高 斯 闭合 法 得 到 的 响应 统计 算 与 
模拟 解 的 结果 吻合 较 好 , 随 着 时 间 演 化 ,位 移 和 速度 
的 方差 逐渐 增 大 并 在 上 = 6s 左 右 到 达 平 稳 。 将 结 
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代入 式 (20) 中 , 取 截 断水 平 n = 6 ,响应 的 概率 密度 

函数 可 以 近似 求 得 ,如 图 3 所 示 。 从 图 可 知 , 不 同时 
刻 的 响应 对 数 概率 密度 函数 近似 解 与 蒙特 卡 洛 模 拟 
解 (MCS) 吻 合 非常 好 ,尤其 在 概率 密度 函数 尾部 。 


这 表明 在 求解 结构 动力 可 靠 性 方面 ,EPC 法 可 以 提 
做 一 个 非常 吊 好 的 概率 密度 函数 近似 解 。 
EPC EMCS 
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四 
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生生 下 人 一 一 MCS 
1 
0 
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| 
加 
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X10 m) 
(b) 4s 
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X/(10 m) 
(c) fF6s 


-0.04 -0.02 0.00 0.02 0.04 


Xi/(10 m) 
(d) 概率 密度 函数 的 演化 曲线 


图 3 平稳 激励 下 不 同时 刻 位 移 的 概率 密度 函数 
Fig.3 PDFs of displacement under stationary 
excitation at different time instants 


位 移 的 概率 密度 函数 演化 曲线 如 图 3(d) 所 示 。 
从 图 可 知 , 随 着 时 间 演 化 ,位 移 概率 密度 函数 曲线 逐 
渐 发 散 。 计 算 结果 表明 在 响应 的 统计 和 矩 达 到 平稳 后 
曲线 的 概率 密度 发 散 程 度 最 大 ,随后 概率 密度 函数 
曲线 保持 不 变 。 

保持 激励 强度 ABR 天] = 0.5 m/s 不 变 , 分 别 
取 平 均 到 达 速 率 和 A = 3,5,20 ,不 同时 刻 的 对 比 结 
如 图 4 所 示 。 从 图 可 知 随 着 到 达 速 率 的 增 大 ,概率 
密度 函数 曲线 逐渐 收缩 。 这 表明 如 果 到 达 速 率 不 是 
足够 大 ,将 该 激励 看 成 高 斯 激励 会 使 得 计算 结果 偏 
于 不 安全 。 


响应 对 数 
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(c) 6s 


图 4 平稳 激励 下 位 移 概 率 密度 函数 与 到 达 速 率 关系 图 


Fig.4 Comparison of PDFs of displacement under stationary 


excitation with different arrival rates 
除 此 之 外 ,不 同 激励 强度 下 结构 的 位 移 概率 密 
度 孔 数 曲线 如 图 5 所 示 。 从 图 可 知 , 随 着 激励 强度 
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的 增 大 ,概率 密度 函数 曲线 逐渐 发 散 。 
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图 5 平稳 激励 下 位 移 概率 密度 函数 与 激励 强度 关系 图 


Fig.5 Relationship between PDFs of displacement under 


-0.06 -0.04 


stationary excitation and different excitation intensities 

现 考 虑 附加 集中 质量 块 对 结构 位 移 的 概率 密度 
函数 的 影响 。 由 表达 式 (7) 可 知 ,改变 和 a 的 值 ， 
仅 引 起 al .as .as 的 改变 ,进而 引起 刚度 y 、 加 速度 
非 线性 参数 8 和 s 的 改变 ,结果 如 图 6 所 示 。 从 图 
可 知 ,附加 质量 主要 影响 结构 的 刚度 y 及 s ,加 速度 
非 线 性 参数 6 变化 范围 较 小 。 因 此 ,便于 分 析 , 取 
wa =1.0 和 oo = 1.0 两 种 情况 研究 附加 质量 对 结构 
响应 的 影响 ,结果 如 图 7 所 示 。 从 图 6 ~7 可 知 , 当 
a =1.0 时 , 随 着 a, 增加 , s 快速 增加 ,位 移 概 率 密 
度 函 数 曲线 呈 发 散 趋势 ; 当 o = 1.0 时 , 随 着 a 增 
加 , s 显著 减 小 , 概率 密度 函数 曲线 出 现 收缩 的 
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0.515 0 
0.680 0 
0.845 0 
1.010 0 
1 17 和 0 
1.340 0 
1.305 0 
1.670 0 
1.835 0 


0.066 84 
0.070 64 
0.071 44 
0.072 23 
0.073 03 
0.073 83 
0.074 63 
0.075 43 
0.076 22 
0.077 02 
0.077 82 


1.720 
2.600 
3.480 
4.360 
5.240 
6.120 
7.000 
7.880 
8.760 
9.640 
10.520 


图 6 结构 参数 与 附加 质量 的 关系 图 


Fig.6 Effects of lumped mass on system parameters 


3.2 ”指数 型 


当 结 构 受 到 非 平稳 泊 松 白 噪声 激励 时 ,通常 激 

励 强度 随时 间 演 化 先 增 大 后 减 小 。 为 了 模拟 这 一 过 

程 ,可 以 采用 NARAYANA 等 "提出 的 指数 型 函数 

作为 时 间 调 制 函数 。 考 虑 一 致 性 ,经 过 Arias 强度 计 

算 , 取 与 平稳 泊 松 白 噪声 激励 持 时 6 s 相同 能 量 的 指 
数 型 函数 

h(t) =18.403(0.0735 +0. 108 3t)exp( —0.056 327) 

(24) 

和 A =3 ,激励 


作为 时 间 调 制 函 数 。 平 均 到 达 速 率 取 
强度 AE[ 产 ] = 0.5 m/s 。 
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首先 ,通过 高 斯 闭合 法 求 得 结构 响应 的 统计 憩 ， 
同样 , 初 值 取 o3(0) = 2o1(0) = 0.001, 如 图 8 所 
示 。 从 图 可 知 ,位 移 和 速度 的 方差 随时 间 演 化 呈现 
ee 并 在 + = 3s 时 达到 最 大 ,随后 
~ a=0.1,a=1.0 
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响应 对 数 
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图 7 位 移 概 率 密度 函数 与 附加 质量 的 关系 图 
Fig.7 Relationship between PDFs of displacement 


and different lumped masses 

利用 统计 矩 结果 ,采用 EPC 法 求 得 位 移 响应 的 
概率 密度 函数 曲线 ,结果 如 图 9 所 示 。 从 图 可 知 , 随 
着 时 间 演 化 ,概率 密度 函数 呈现 先 发 散 ,后 收缩 的 趋 
势 , 并 且 在 上 =3s 时 发 散 程 度 最 大 ,这 与 响应 的 方差 
演化 趋势 相同 。 

考虑 不 同 的 到 达 速 率 对 结构 响应 的 影响 ,分 别 
取 A =3,5,20 ,保持 激励 强度 不 变 , 以 1 = 3 s 为 例 ， 
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位 移 啊 应 的 概率 密度 函数 如 图 10 所 示 。 从 图 可 知 ， 
es 啊 应 计算 结果 相同 , 随 着 到 


速率 的 增 大 ,响应 的 概率 密度 函数 曲线 逐渐 收缩 。 
QO E[x’ ](EPC) 
0.005 ; xz](EPC) 
—-E[x:](EPC) 
0.004 ——E[x’](MCS) 
际 — Elxx,](MCS) 
0 Elx:](MCS) 
NY 
0.002 
0.001 bi 
0.000 
0 40 80 120 160 200 
1W(10 s) 
图 8 非 平 稳 激 励 下 结构 响应 的 统计 算 
Fig.8 Moments of system responses under 


nonstationary excitation 


响应 对 数 
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图 9 非 平稳 激励 下 不 同时 刻 位 移 的 概率 密度 函数 
Fig.9 Probability density functions of displacement under 


nonstationary excitation at different time instants 


-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 
X/(10 m) 


图 10” 非 平稳 激励 下 位 移 概 率 密度 函数 
与 到 达 速 率 关 系 图 
Fig. 10 Relationship between PDFs of displacement under 
nonstationary excitation and different arrival rates 

最 后 ,考虑 不 同 激励 强度 对 结构 响应 的 影响 , 结 
果 如 图 11 所 示 。 从 图 可 知 ,与 平稳 激励 下 的 情况 相 
同 , 随 着 激励 强度 的 增加 ,结构 位 移 的 概率 密度 函数 
呈 发 散 趋势 。 


微 信 公 众 号 :应 用 力学 学 报 


900 应 用 力学 学 报 


3. 
| AE(Y’)=0.1 m/s “AE(Y 0.3 mY/s 
-AE(Y )0.5 m/s AE(Y =0.7 m/s 
1— 180°1.0 m/s 
渡 0 
长 
| -1 
- 
= 多 
3 
-未 有 | : oa, 
-0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 0.08 
X10 m) 


图 11 非 平稳 激励 下 位 移 概率 密度 函数 与 激励 强度 关系 图 
Fig. 11 Relationship between PDFs of displacement under 


nonstationary excitation and different excitation intensities 


4 结 论 


本 研究 首先 根据 Lagrange 方程 建立 了 非 高 斯 激 
励 下 附加 集中 质量 块 的 悬臂 梁 随 机 微分 方程 ,通过 
引 数 多 项 式 闭 合法 求解 了 在 平稳 泊 松 日 噪声 激励 和 
非 平 稳 泊 松 白 噪声 激励 下 结构 响应 的 概率 密度 函 
数 ,并 分 析 了 平均 到 达 速 率 ` 激 励 强度 和 附加 集中 质 
量 对 结构 的 影响 ,主要 结论 有 :1) 采 用 指数 多 项 式 闭 
合法 计算 的 结果 与 模拟 解 吻合 较 好 ,尤其 是 对 可 靠 
性 分 析 起 重要 作用 的 尾部 ;2 ) 平 均 到 达 速 率 增 大 , 结 
构 响 应 的 概率 密度 函数 曲线 逐渐 收缩 ,因此 ,将 外 部 
激励 简单 看 作 高 斯 激励 会 使 得 计算 结果 侦 于 不 安 
全 ;3) 激 励 强 度 越 大 ,结构 响应 的 概率 密度 函数 越发 
散 ;4) 附 加 集中 质量 块 的 质量 保持 不 变 , 增 大 转动 惯 
量 , 结 构 响 应 的 概率 密度 函数 曲线 快速 收缩 ;保持 转 
动 惯量 不 变 , 增 大 质量 ,结构 响应 的 概率 密度 函数 曲 
线 呈 发 散 趋势 。 
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